Conservation laws for the voter model in complex networks 
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We consider the voter model dynamics in random networks with an arbitrary distribution of the 
degree of the nodes. We find that for the usual node-update dynamics the average magnetization 
is not conserved, while an average magnetization weighted by the degree of the node is conserved. 
However, for a link-update dynamics the average magnetization is still conserved. For the particular 
case of a Barabasi- Albert scale-free network the voter model dynamics leads to a partially ordered 
metastable state with a finite size survival time. This characteristic time scales linearly with system 
size only when the updating rule respects the conservation law of the average magnetization. This 
scaling identifies a universal or generic property of the voter model dynamics associated with the 
conservation law of the magnetization. 

PACS numbers; 64.60.Cn,89.75.-k,87.23.Ge 



I. INTRODUCTION 

Conservation laws play an important role in the char- 
acterization and classification of different nonequilibrium 
processes of ordering dynamics. For example, in Ki- 
netic Ising models one distinguishes between Glauber and 
Kawasaki dynamics. In Glauber dynamics the individ- 
ual dynamical step is that of flipping a spin, while in 
Kawasaki dynamics two nearest neighbor spins are ex- 
changed. Kawasaki dynamics conserves magnetization 
and Glauber dynamics does not. As a consequence, the 
Glauber and Kawasaki dynamics give rise to different 
scaling laws for domain growth in coarsening processes 
0, and they define different nonequilibrium universality 
classes. The conservation law of Kawasaki dynamics is 
implemented at each time step of a stochastic dynamics. 
A different type of conservation law is the one that refers 
to an ensemble average. An example of such conserva- 
tion laws is the conservation of the average spin (global 
magnetization) in the voter model HH. When studying 
spin dynamical models in regular lattices the existence of 
an ensemble conservation law does not imply an elemen- 
tary step conservation such as imposed in the Kawasaki 
dynamics. Recent interest in ordering processes focuses 
in situations in which the spins are located in the nodes 
of a complex network, i.e. a network with a large het- 
erogeneity in the number of nearest neighbors with which 
each spin interacts. This does not affect the fulfillment of 
a conservation law of the type of the Kawasaki dynamics, 
but the implementation of an ensemble average conserva- 
tion law requires a careful thought of the dynamical rules. 
As an interesting example, we discuss in this paper the 
differences between node and link-update dynamics for a 
voter model in which spins are located in the nodes of 
a random network with an arbitrary degree distribution. 
Only link-update dynamics respects the global magneti- 
zation conservation law, while another conservation law 
exists for node-update dynamics. 

The standard voter model is defined by a set of "vot- 



ers" with two opinions (spins Ui ~ ±1) located in the 
nodes of an hypercubic lattice. The elementary dynam- 
ical step consists in randomly choosing one node (asyn- 
chronous update) and assigning to it the opinion of one 
of its nearest neighbors, also chosen at random (node- 
update). One time step corresponds to updating a num- 
ber of nodes equal to the system size, so that each node 
is on average updated once. In d = 1 this dynamics is 
equivalent to the zero temperature Glauber kinetic Ising 
model. In general dimensionality, the global magneti- 
zation is conserved in the thermodynamic limit of large 
systems and the dynamics is dominated by interfacial 
noise. The infinite system coarsens for d < 2, with a 
slow logarithmic decay in the critical dimension d — 2. 
At variance with other ordering dynamics, coarsening 
takes place here without surface tension. The role of 
the conservation law of the magnetization and of the Z2 
symmetry (±1 states) in the voter dynamics universality 
class has been studied in detail in the critical dimension 
d = 2 of regular lattices j^. We are here interested in 
situations in which there is no long-time coarsening: for 
d > 2 a, finite system in a regular lattice reaches one of 
the homogeneous attractors in which all the A'' spins have 
the same value. The time to reach such consensus t, or 
survival time, scales as r oc iV, so that there is no com- 
plete ordering in the thermodynamic limit . This same 
scaling behavior has been also found for the voter model 
in a small- world network 0, 13 so that it can be identi- 
fied as a generic property of the voter model dynamics. 
We show in this paper that for a Barabasi- Albert (BA) 
scale-free network Q the scaling law r oc is only ob- 
tained when the average magnetization is conserved, that 
is when link-update dynamics is used. 



II. NODE VS LINK-UPDATE IN THE VOTER 
MODEL 

Generically, in a complex network such as the small- 
world network, there is an heterogeneous degree distri- 
bution with nodes having a different number of links. In 
this case, as explained in more detail below, node-update 
dynamics does not guarantee the conservation of the av- 
erage magnetization. The conservation is guaranteed in 
a link-update dynamics in which the elementary dynam- 
ical step consists in randomly choosing a pair of nearest 
neighbor spins, i.e. a link, and randomly assigning to 
both nearest neighbor spins the same opinion when they 
had different opinion, and leaving them unchanged oth- 
erwise. The reported simulations of the voter model in a 
small world network seem to use a node- update , while 
the analytical results are obtained in an approximation 
which enforces conservation of the average magnetiza- 
tion 8]. Since in both cases the scaling law r cx iV is 
obtained, the role of the conservation law in this generic 
property is unclear. Still, the degree heterogeneity in a 
small world network is rather small, and therefore its ef- 
fect on the conservation law for node-update is probably 
not significant. The question is much more crucial when 
considering, for example, a scale-free network which has 
nodes with large degree heterogeneity, as we do next. 

III. CONSERVATION LAWS 

Results for the ensemble average normalized magneti- 
zation, 

{<^it)) = ^ — — ^ , (1) 

of a voter model dynamics in a BA network Q are shown 
in Fig.nfor node-update and link-update dynamics. The 
ensemble average indicated as (•) is an average over the 
realizations of the stochastic dynamics and different ini- 
tial conditions. For both updating rules we have consid- 
ered the same initial distribution of spins in the initial 
configuration, with half of the nodes with spin +1, and 
half with spin —1. To provide more clear evidence of the 
dynamical differences of the two update rules we have 
chosen an initial configuration in which all the nodes that 
were given initial spin -f 1 have a higher degree than those 
that were given initial spin —1. For the node- update dy- 
namics we find that the average spin is not conserved, 
but rapidly changes towards the positive side, due to 
the influence of the high-degree nodes in their immediate 
neighborhood. For the link-update dynamics the average 
magnetization remains at its vanishing initial value. 

For any complex network with an heterogeneous dis- 
tribution of the degree of the nodes, the differences be- 
tween the two updating rules are easily understood as 
follows. The conservation of the average magnetization 
in a regular lattice relies on the fact that, if node i and 
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FIG. 1: Ensemble average magnetization (a) in a BA net- 
work. The most connected half of the nodes in the network 
have initial spin ai = -1-1, while the other half has initial spin 
CTi — —1. System size N = 1000, average degree k — 8, and 
averaged over 1000 realizations. Circles: node-update; Dia- 
monds: Link-update. 



j have different opinions and are connected, the proba- 
bility for j to change to the spin of i in a time step of 
the dynamics is the same as for i to change to the spin 
of j. Since this is true for all nodes, the ensemble aver- 
age magnetization is conserved. However, in a network 
where the nodes have not the same degree, this is not 
longer true for node-update dynamics. For instance, if 
« is a highly connected node with degree ki, and j is a 
node with a low degree kj < ki, and they have different 
spins, then the probability Pij that i changes to the spin 
of j {{Nki)"-^) is smaller than the probability Pji that 
j changes to the spin of i {(Nkj)~^). Thus, the average 
magnetization is not conserved. Choosing the node to 
be updated preferentially, so that Pji = Pij , makes the 
average spin conserved again. Preferentially choosing the 
node to be update in this way is equivalent to randomly 
choosing a link in the network and updating it in random 
direction (link-update) . 

As we have argued the ensemble average normalized 
magnetization (a) is not conserved in a complex network 
in the voter model with node-update rule. In order to 
compensate for the different degrees of the nodes, we con- 
sider a degree weighted normalized magnetization [Tol | 

where ki is the degree of the node i. The total number of 
links is introduced in the denominator to normalize the 
weighted magnetization between [—1, 1]. 

If we define S{t) = 'Y^^^-^^kiai{k,t), then for a given 
configuration and averaging over stochastic realizations, 
the expected change (ASij) in a time step of the dy- 
namics due to node i changing spin to the spin of node j 




FIG. 2: Time evolution of {a{k,t)) for node-update dynam- 
ics in a BA network with k — 8. Initial configuration as in 
Fig. System size A'^ = 10, 000, average degree k — 8 and 
average taken over 10, 000 realizations. Thin lines correspond 
to numerical data; thick lines correspond to the analytical 
predictions (Eq. 110|l 'l for the initial values a{k,0) = —1 and 
+ 1. 



is given by 



while 
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where (•)^ represents an average over stochastic real- 
izations for a given configuration c at time t. Thus 
(AS'ij -f AS'ji)^ — 0. As this argument apphes for any 
pair of neighbors, and any configuration, the ensemble 
average weighted magnetization (S) is conserved for the 
voter model using node-update in complex networks with 
arbitrary degree distributions. 

Given the conservation law of (S) we can find the 
asymptotic value of the average magnetization. To this 
end, we introduce the normalized magnetization of the 
nodes with given degree k as: 



a{k,t) 



(5) 



where is the number of nodes with degree k and the 
sum in the numerator is over all nodes with the same 
degree k. 

For a given configuration c, the expected change of 
spin of a node i with degree ki due to interaction with 
its neighbors in a time step of the dynamics is given by 
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(T.j{t) - a^(t) 



(6) 



10^ 



6 8 10 12 14 
average k , , 



10^ 



N 



10 



10 



FIG. 3: Survival time for node-update dynamics in BA net- 
works of different sizes A'', and average degree fc = 6. Inset: 
Survival times for different average degree. All data obtained 
from at least 1000 realizations of networks of size 10000. 



where Vi is the neighborhood of node i, that is the nodes 
connected by a link to node i. From this expression, if 
we add for all the nodes with the same degree we obtain 



{Aa{k,t))= J2 



i:ki=k jeVi 



a J (t) - cTj (t) 



(7) 



We can now split the r.h.s of Eq. ([TJl in two terms. The 
second term is simply cr(fc,f). For the first term, we as- 
sume a mean field approximation, i.e., we consider a ran- 
dom network where the sum over neighbors is equivalent 
to a random sampling over the whole network. Then 



i:k,=kjeV^ 



EkPik)kaik,t) _ l:^=l,Nk^O^{t) 
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J2i=l,N 



(8) 

where P{k) = Nk/N is the degree distribution of the 
network, that is, the probability to find a node of degree 
k. 

Thus, (A(T(fc,t)}^ = T,-a{k,t), where T.{t) and cr(fc,i) 
are calculated in the given configuration c. Averaging 
over different configurations, we find the evolution equa- 
tion of the ensemble average oia{k,t) [13: 



d{a{k,t)) 
dt 



{a{k,t)) 



(9) 



Given the conservation law for (E), the solution of 
Eq. (PI is an exponential approach to to the asymptotic 
value 



(a(fc,i)) = ((a(fc,0))-(E))e-* + (E) 



(10) 



In the long time limit {a{k,t)) approaches a constant 
value (E) independent of k. Therefore this constant value 
coincides with the long time limit of the ensemble average 
normalized magnetization {a{t — > oo)). This value is a 
property of the ensemble of initial configurations. 



A numerical check of these general results for the par- 
ticular case of the BA network is shown in Fig. [21 where 
the fast exponential decay to the final average value is 
shown. The asymptotic value corresponds to the analyt- 
ical prediction of (E) calculated in the initial configura- 
tions, that is, the ensemble of initial configurations. 

IV. SURVIVAL TIMES IN BARABASI-ALBERT 
SCALE-FREE NETWORKS 

We next study the consequences of the two different 
updating rules and associated conservation laws in the 
ordering dynamics of the voter model in a BA network. 
For any of the two updating rules the system falls, af- 
ter an initial transient, in a metastable partially ordered 
state until a finite system size fluctuation takes the sys- 
tem to one of the two ordered attractors of the dynamics. 
This qualitative behavior is similar to the one found in a 
small world network and like wise can be characterized 
in terms of the temporal evolution of the average inter- 
face density p, defined as the density of links connecting 
sites with different opinions. In a given realization of the 
dynamics p initially decreases indicating a partial order- 
ing of the system. After this initial transient p fluctuates 
randomly around an average value until a fluctuation or- 
ders the system leading to an absorbing state with p — Q. 
Considering an ensemble of realizations, the ordering of 
each of them happens randomly with a constant rate. 
This is reflected in an exponential decay of the average 
interface density 

{p)<xe-i , (11) 

where r is the survival time of the partially ordered 
metastable state. This survival time turns out to be a 
quantity that diverges with growing system size TV, so 
that the system does not order in the thermodynamic 
limit. 

We have measured the characteristic time r for the 
two updating dynamical rules and for different system 
sizes N and different mean degree k. Our results are 
summarized in Figs. Eland 0] For the node-update rule 
(Fig. O, in which there is no conservation law of the 
average magnetization, we have found that r scales with 
system size N as 

T oc , (12) 

where 7 — 0.88 ± 0.01. For a fixed system size, the value 
of the average degree k does not seem to have any def- 
inite influence on r (see the inset of Fig. O. The value 
of the exponent 7 is consistently and significantly differ- 
ent from 7 = 1 which is the exponent analytically found 
for regular hypercubic lattices and for an annealed small 
world network (8| . In these last two cases the dynamical 
rules respect the conservation of the average magnetiza- 
tion. When we implement the link-update rule which 
also conserves the average magnetization we find a result 
for 7 which is consistent with 7=1 (Fig. 
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FIG. 4: Survival time for link-update dynamics in BA net- 
works of different sizes N , and average degree k = 8. AH data 
obtained from 1000 realizations. 



V. CONCLUSIONS 

In summary, we have shown that the voter model dy- 
namics does not lead to an ordered state in a scale-free 
network in the thermodynamic limit. This is consistent 
with the results for a small world network, and in general 
for networks of dimensionality d > 2. Finite size effects 
order the system in a time which depends on the updat- 
ing dynamical rule. Only for the updating rule fulfilling a 
conservation of the global magnetization this time scales 
linearly with the system size. This is also consistent with 
the result for regular hypercubic lattices of d > 2 and for 
the voter model in annealed small- world networks 0. 
Such scaling can then be taken as a proper characteriza- 
tion of universal properties of the dynamics of the voter 
model. 

We note that there are several instances in which the 
conservation law of the global magnetization is naturally 
broken independently of the updating rule, as for exam- 
ple the consideration of a zealot "12] , or the dynamics in 
a directed network. On the other hand there are spread- 
ing phenomena with "spins" having N > 2 states which 
show voter-like generic properties in d=2 regular lattices 
and that fulfill the global magnetization conservation law 
[T^ . Further study of these other cases and the impli- 
cations of other conservation laws, as the one reported 
here for node-update dynamics, would be useful to iden- 
tify generic and nongeneric properties of the voter model 
dynamics. 
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